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Abstract 

The role of BEC for photons and neutral pions is briefly discussed and 
illustrated by means of new method of numerical modelling of BEC proposed 
by us recently. 
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Photons and vr^'s play particular role in high energy multiparticle processes, es- 
pecially in heavy ion collisions, which are believed to be the place where the Quark 
Gluon Plasma is possibly formed. This is because they are immune to the Coulomb 
final state interactions distorting BEC analysis of charged particles. The vr'^ inter- 
ferometry can be then used as a kind of check allowing to estimate the correcteness 
of procedures accounting for this effect in BEC between charged particles. Of minor 
experimental consequensce at present (but otherwise very interesting) is fact that, 
because both photons and vr^'s are also their own antiparticles, one expects that for 
low momenta (actually, for k^^^^o —>■ 0) one should observe some additional effects 
emerging from the particle-antiparticle correlation, a pure quantum mechanical ef- 
fect discussed at length in jTI. Photons are also the only particles which could, in 
principle, provide us with detailed picture of the whole history of heavy ion col- 
lision process because they are produced anywhere where one has deceleration or 
acceleration of charges whereas hadrons (including vr^'s) are produced only at final 
hadronization period of collision. This means that they are produced starting with 
very early stage of collision (sometimes called " preequilibrium" ) , through the stage 
of formation and equilibration of the QGP and through the mixed stage where first 
hadrons are being produced in QGP environment and ending with the last kinetic 
freez-out stage where hadrons are finally formed (see |2j for most recent review; no 
BEC discussion is presented there, however). It is obvious that production mech- 
anism of photons at different periods of collision will be different, ranging from 
bremstrahlung at the beginning, via thermal emission at QGP phase, to decays of 
neutral particles (mostly 7r°'s) at the end For us the most interesting would be 
directly produced photons, especially thermal ones originated from the QGP. How- 
ever, according to the recent WA98 data [3], they consist of only a few percent of 
the all photons visible. 

Photons should be also copiously produced at PHOS detector of ALICE at LHC, 



BEC measurements, both with direct and decay photons (for references see 0, the 
most recent paper on this subject is j3] where updated references can be found). 
Whereas the method to be used for BEC are essentially the same as for other parti- 
cles (cf. IHC]) there are some new problems one is facing which are connected with 
the fact that at ALICE one probably will not be able to measure directly tt'^'s but 
only photons they are decaying to. Question which arises then is: notwithstand- 
ing this fact, can one deduce BEC characteristics of vr^'s from BEC characteristics 
of 7's? This problem has been addressed in |6j and positive conclusion has been 
reached as well as the procedure of such analysis has been provided, which can eas- 
ily be applied when data will be available. However, opposite conclusion has been 
reached in j3] where photons from vr^'s were demonstrated to show very character- 
istic structure at small values of Q (of the order of pion mass) and fiat distribution 
for larger values of Q. On the other hand, it was also shown there that, at least 
for SPS and RHIC energies, one can actually separate correlations of photons pro- 
duced directly in the hot zone and residual correlations of decay photons (i.e., those 
coming predominantly from tt^'s). It turns out that (at least in the framework of 
the hydrodynamical model of collision used there) the former contribute essentially 
to the region of small invariant relative momenta, Qinv ^ 50 MeV, while the later 
dominate the large relative momomentu. 

Here we shall present numerical calculations in which vr^'s were created by Monte 
Carlo algorithm, which incorporates already BEC effect Q\. Referring to |7j for 
details and physical motivation we shall concentrate here on the same problem as 
was considered analytically in [Sj and jB], namely, to what extent BEC picture of 
original vr^ is visible in finally observed 7's to which they decay and what is the 
chance to detect BEC of originally produced 7's as well. To this end we have to 
modify slightly algorithm used in j2j by allowing for varying transverse momenta 
of produced vr^'s (they were fixed as px = {vt) in fZl, i-e., the case considered was 



choosing energy E of selected particles (done in the same way as before j7j) we 



shall notice that E = yp^ + pi + /i^ (where n denotes pion mass) and choose now 
Pt = IptI from distribution P{pt) = gw{~Pt / {pt)) (in calculations presented here 
parameter (px) = 0.2 GeV/c) such that pr € (0., \/E^ — /i^). Then the azimuthal 
angle G [0,27r] was selected from uniform distribution and direction "right-left" 



for longitudinal momentum pi^ was chosen. In this way pi^ = ±y E'^ ~ Pt ~ /^^- 
this way we star a new elementary emitting cell (cf. |0] for more details) and add to 
it other particles of the same kind momenta of which were selected from gaussian 
distributions with widths and ax treated as free parameters. It turnes out that in 
the example considered here results are not particularly sensitive to their difference 
therefore we have put them equal: ai = (^t ^ ■ 
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Fig.l Example of BEC results for tt" as produced by our algorithm (full lines) and BEC for 7's resulting from their 
decays. Left panel shows results for different sizes of elementary emitting cells (EEC's) characterized by parameter 
CL.T- Right panel shows results for different strength of BEC effect characterized by parameter Pq (leading to 

different occupancy of the average EEC, cf. 7 ). 

Our results are shown in Fig. 1. They confirm results of 0, namely that photons 



^The reason is that our example is still very much one-dimensional. The reason is that in |7j we 
wanted to check whether it is possible in principle to reproduce e + e~ annihilation data and their 
phase space is clearly a quasi-one-dimensional. However, this limitation does not affect results and 



from decaying pions do not follow their BEC pattern but show, instead, their own 
pattern, which is quite robust against different choices of parameters, and which 
consist in essentially constant (but small) correlations for small relative momenta 
vanishing after Q of the order of pion mass. The possible explanation of this result 
seems to be the same as in (modulo possible deviations in numbers caused by 
our limitation in px and essentially 3-dimensional case considered in 
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^At this moment we cannot offer any explanation of difference between results of |S1 and [Sj 
noticing only that we do not think that explanation provided in |0| and attributing this difference to 
different way of introducing cuts in momenta should not work for the first panel of figure presented 



